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Abstract: We continue our study of the correlators of a recently discovered family 
of BPS Wilson loops in M = 4 supersymmetric U(N) Yang-Mills theory. We per- 
form explicit computations at weak coupling by means of analytical and numerical 
methods finding agreement with the exact formula derived from localization. In par- 
ticular we check the localization prediction at order g 6 for different BPS "latitude" 
configurations, the Af = 4 perturbative expansion reproducing the expected results 
within a relative error of 10~ 4 . On the strong coupling side we present a supergrav- 
ity evaluation of the 1/8 BPS correlator in the limit of large separation, taking into 
account the exchange of all relevant modes between the string worldsheets. While re- 
producing the correct geometrical dependence, we find that the associated coefficient 
does not match the localization result. 
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1. Introduction 

The maximally supersymmetric Yang-Mills theory provides the simplest dynamics 
among the four- dimensional gauge theories and has represented an important and 
interesting laboratory from the theoretical perspective. Not only is it believed that 
the theory has an exact dual description, the type IIB ten-dimensional string theory 
in the AdSs x S 5 background |l[], but new fascinating connections have appeared 
throughout the years. The integrable structures underlying the spectrum of the 
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anomalous dimensions 0, |3| , the exact exponentiation properties observed in scatter- 
ing amplitudes Q and the possible connection with the geometric Langlands program 
[[J are just a few examples of the richness still hidden in M = 4 SYM. 

Remarkably exact results exist also for Wilson loops, which in M = 4 theory can 
be generalized to preserve some amount of superconformal symmetry. The simplest 
operator of this kind is a circular Wilson loop which couples to one of the six adjoint 
scalar fields of the theory: it is called the 1/2 BPS Wilson loop because it preserves 
one half of the 32 superconformal symmetries. It was conjectured || [7| that the 
expectation value of such an operator can be computed in the Gaussian matrix model. 
The conjecture was supported by an explicit two-loop perturbative computation, 
while from the dual string theory point of view, in a suitable limit of large N and 
large 't Hooft constant A = gy-M^i ^ ne Gaussian matrix model nicely agrees with the 
solution to the minimal area problem ||, |9|. More generally other kinds of Wilson 
loops, which preserve various amounts of supersymmetry have been constructed and 
studied. In particular a family of 1/16 BPS Wilson loops of arbitrary shape on 
a three-sphere S 3 , embedded in the Euclidean four dimensional space-time, were 
presented in |L0|, |TTJ . Restricting the contour of the loops to the equator one gets 1/8 
BPS Wilson loops and a conjecture in this case has been proposed: the expectation 
value of such Wilson loops is captured by the zero-instanton sector of the ordinary 
bosonic two-dimensional Yang-Mills living on the S 2 . The coupling constant of the 
2d Yang-Mills theory is related to the coupling constant of the A/"=4 SYM as g\ d = 
—gl d /(27rr 2 ), where r is the radius of the S 2 . The loops are still computed by a 
Gaussian matrix model, the circular Wilson loop being a particular case of this 
general family. 

This conjecture was further supported at order g 4 in fl2], [L3|] for the expectation 
value of a single Wilson loop operator of arbitrary shape on S 2 , and further extended 
at the level of BPS correlators of Wilson loops |i~4], |i"5fl, where the multi-matrix model 
describing the correlators in the zero-instanton sector of YM2 has been derived. 
Non-trivial consistency checks and explicit computations supporting the conjecture 
for correlators have also been performed ||15|| . The conjecture has been recently 



extended to include 't Hooft loops and S-duality, taking into account non-trivial 
instanton sectors [IB|. 



The emergence of a two-dimensional theory underlying the dynamics of some 
BPS sectors of Af=4 SYM can be understood from the localization of the four- 
dimensional path integral to particular supersymmetric configurations, as recently 
shown by [[17], [Tj|. The moduli space of solutions to the supersymmetry equations 
is parameterized by two-dimensional data and the effective action governing the rel- 
evant dynamics is the semi-topological Hitchin/Higgs- Yang-Mills theory: the com- 
putation of regular 1/8 BPS observables can be mapped there and reduces to usual 
YM 2 on S 2 . 

The localization of the path-integral in four-dimensional supersymmetric gauge 
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theories is not a novelty: the exact computation of the prepotential in M = 2 SYM 
has been derived in |19| through this kind of procedure, summing up all instan- 
ton contributions. Actually the result concerning the 1/2 BPS Wilson loop can be 
extended to TV = 2* SYM theories, taking also into account the contribution of 
instantons (which decouple in the M = 4 case). Quite recently, from the exact ex- 
pression of the partition function derived from the localization approach of [I7|, a 
rather general class of Af = 2 superconformal gauge theories introduced in has 
been shown to be described by two-dimensional Liouville theory f21 |. 

It appears therefore important to test the results expected from path-integral 
localization through the familiar perturbative QFT methods and the AdS / CFT cor- 
respondence, even in the simplest case of M = 4 where some checks are still missing. 
In particular the two-dimensional gauge theory should compute not only the ex- 
pectation value of a single 1/8 BPS Wilson operator but even correlators of loops 
preserving the same amount of supersymmetry. A first step in this direction was 
taken in [Tj|, where an apparent disagreement was observed in the limit of coinci- 
dent loops. Later the relevant matrix-model result was shown to be consistent with 
the supergravity picture [Q. In [15 we started a systematic approach to the com- 
putation of the correlator of two "latitude" BPS Wilson loops, at weak coupling by 
perturbation theory and at the strong coupling through AdS/CFT correspondence. 
We checked the formula derived from the zero-instanton sector of YM2 at order 
0(g A ) and we showed that, in the limit where one of the loops shrinks to a point, 
logarithmic corrections in the shrinking radius are absent at 0(g 6 ). This last result 
strongly supported the validity of the general expression and suggested the existence 
of a peculiar protected local operator arising in the OPE of the Wilson loop (see also 
^] for a related investigation). Using the string dual of the M = 4 SYM correlator 
in the limit of large separation, we also presented some preliminary evidence for the 
agreement at strong coupling. 

In this paper we continue our study of the two-latitude correlator, extending 
our previous investigations. First of all we present strong evidence that the weak 
coupling perturbative computation agrees with the matrix model expression. We 
evaluate numerically the expectation value of correlators at order g 6 for two particular 
configurations: a symmetric one in which the loops are two latitudes at polar angle 
= 5 and 9 = 7r — 5 (1/4 BPS system) and the other with a loop fixed on the 
equator 6 = n/2 and the second at generic angle 9 = 5 (1/8 BPS system). No 
particular limit has been considered and the agreement is quite good over the whole 
range of our study, including angles 5 between 0.7 and tt/2. The relative error 
between the YM 2 prediction and the f\f = 4 SYM calculation is of the order of 10~ 5 . 
For values of 5 less than 0.7 the requirement on the precision of the calculation of 
certain integrals becomes prohibitive. Generically, the errors involved grow in the 
opposite (coincident) limit of 5 = tt/2, however we find that they are manageable 
even when 5 is very close to tt/2. The supergravity calculation is also tackled and 
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should reproduce the strong coupling result, at large N, of the exact localization 
answer: unfortunately we were not able to find such agreement. We compute the 
exchange of supergravity modes between the widely separated worldsheets describing 
the Wilson loops at strong coupling. We identify all modes contributing to the 
correlator at leading order in the large separation limit. While the sum of these 
exchanges produces a qualitative agreement with the matrix model, we observe a 
deviation in the numerical coefficient. We comment on this puzzle and will leave its 
resolution to future investigations. 

The plan of the paper is the following: in Section 2 we briefly recall the structure 
of the BPS Wilson loops, their expectation values and correlators according to the 
localization formula and discuss our previous results. In Section 3 we present our 
numerical computation in detail, explaining our procedure and critically examining 
the numerical agreement. In Section ^ the strong coupling computation is performed 
in detail using the familiar methods of AdS/CFT, and the origin of the mismatch 
is discussed. In Section 5 we draw our conclusions and discuss future directions of 
research. 



2. The supersymmetric Wilson loops and their correlators 

We start by considering the family of BPS Wilson loops that has been introduced in 
[LlJ : a simple way to understand this construction is to observe that it is possible to 



pack three of the six real scalars present in TV = 4 SYM into a self-dual tensor 
and to use the modified connection 

— > An + i<S>„ v x v (2.2) 

in the Wilson loop. The crucial elements in this definition are the tensors c* : they 
can be defined by the decomposition of the Lorentz generators in the anti-chiral 
spinor representation (7^) into Pauli matrices 7$ 

\{l-l 5 )l» v = i°U T ii ( 2 - 3 ) 

where the projector on the anti-chiral representation is included (7 s = — 7 1 7 2 7 3 7 4 )- 
The matrix M l j appearing in ( |2.1| ) is 3 x 6 dimensional and is norm preserving, i.e. 
MM T is the 3 x 3 unit matrix (an explicit choice of M is M 1 i = M 2 2 = M 3 3 = 1 
and all other entries zero). 

More geometrically, the tensors <t^„ are related to invariant one-forms on S 3 

a R,L = 2 [ ± ( x 2 dx 3 _ x 3 dx 2^ + ^A dx l _ x l dx ^ 

a*' L = 2 [±{x 3 dx l - x l dx 3 ) + {x 4 dx 2 - x 2 dx 4 )] (2.4) 
2 [±(x 1 dx 2 - x 2 dx 1 ) + (x 4 dx 3 - x 3 dx A )] , 



- 4 - 



where erf are the right (or left-invariant) one- forms and erf are the left (or right- 
invariant) one-forms: explicitly 



erf = 2o i lxv x»dx v . 



;2.5) 



The BPS Wilson loops can then be written in terms of the modified connection 



Aft + i$ uu x u as 



W = -jy Tr V exp I dx» (iA^ - a\ v x v M i /$ J ) . 



;2.6) 



Actually the operator 



is supersymmetric only when the loop is restricted to a 



(2.7) 



three dimensional sphere. This sphere can be taken to be embedded in M , or as a 
fixed-time slice of S 3 x R. The authors of [JTl ] have shown that requiring that the 
supersymmetry variation of these loops vanishes for arbitrary curves on S 3 leads to 
the two equations 

7Mi + ^P*7 5 eo = , 

7Mq + ^A 5 ei = , 

that can be solved consistently: for a generic curve on S 3 the Wilson loop preserves 
1/16 of the original supersymmetries. We remark that this construction needs the 
introduction of a length-scale, as seen by the fact that the tensor ( |2.1| ) has mass 
dimension one instead of two: we will fix the scale to be the radius r of S 3 . 

The situation becomes more interesting for special curves, when there are extra 
relations between the coordinates and their derivatives: in this case there will be 
more solutions of (\2.7j) and the Wilson loops will preserve more supersymmetry. A 
particularly interesting case is when the loop lies entirely on a S 2 : it is possible to 
show that these operators are generically 1/8 BPS and Wilson loops lying on the 
same two-sphere enjoy common supersymmetries. Inspired by the explicit evaluation 
of the first non-trivial perturbative contribution the authors of |Tl[ conjectured that 
the 1/8 BPS Wilson loops constructed on S 2 can be exactly calculated, claiming the 
equivalence with the computation of Wilson loops in ordinary YM 2 on the sphere, in 
the zero-instanton sector [23]. Yang-Mills theory on a Riemann surface is completely 
solvable |24|] and the exact expression for the Wilson loop is also available [^5|: the 



restriction of the full answer to the zero-instanton sector follows from rewriting the 



exact solution as an instanton expansion [26]. Based on this relation, the following 
exact formula for the quantum expectation value of the 1/8 BPS Wilson operator 



(W) 



N 



_ 2 A 1 A 2 
94d A 2 



exp 



2 A 2 



:2.8i 



was proposed in [11|, where Lj v „ 1 (x) is a Laguerre polynomial, A is the area of the 



sphere and A\^ are the areas enclosed by the loop. The result follows by identifying 
the two-dimensional coupling constant g\ d with the four-dimensional one through 
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g\ d = —g\ d /A and is equal to the expectation value of the circular Wilson loop, 
which is computed in a gaussian Hermitian matrix model |], [7|, 



(W c ) 



N 



Trexp(M) 



1 

Z 



VM—Ty 

N 



exp(M) 



exp 



2~TrM z 



(2.9) 



after a rescaling of the coupling constant g\ d — * g\ d A\A 2 j A 2 '. The conjecture was 
further supported at the second non-trivial perturbative order gf d in [12, 13 1 for the 



expectation value of 1/8 BPS Wilson loop operators of various shape on S 2 while the 
emergence of the two-dimensional theory underlying this peculiar dynamics in A/"=4 
SYM has been understood from the localization of the four-dimensional path integral 



to particular supersymmetric configurations |L7|, |lq| - According to this procedure, 
the computation of M = 4 observables through Yang-Mills theory on S 2 depends 
just on the presence of some preserved supersymmetry: correlators of 1/8 BPS loops 
lying on the same sphere should therefore be computable as well in terms of the zero- 
instanton sector of two-dimensional Yang-Mills. The relevant correlators have been 



derived in [|14], [L5J and are easily obtained from a multi-matrix model: disregarding 
instanton contributions, the formula for the correlator of two BPS loops winding 
respectively n\ and n 2 times around themselves is 



W{A 1 ,A 2 ) 



C N N 2 



DV 1 DV 2 e 



2 92d A T- A 3 



_A^+A 2 _ 

2 !>'L A 2 A 3 



Tr(V 2 2 )+ 



s'L A 3 



Ti(e iniVl )Ti(e m2V2 ) 



Tr(yiy 2 ) 



in2V2\ 



(2.10) 



where the normalization is chosen to be 



C 



N 



DVxDV 2 e 2g ™ 



Al + A 3 , 2) A 2 +A 3 T /y2u 1 ^yy) 



2 4d A l A 3 



(2.11) 



and A 3 = A — A\ — A 2 . The very same result has also been obtained from Feyn- 
man graph calculations using the Mandelstam-Leibbrandt prescription for the vector 
propagator in light-cone coordinates and resumming perturbation theory to all or- 
ders p2| . The final matrix integrals can be exactly computed at finite iV in terms of 
Laguerre polynomials JT5] . For small g 2 d this expression can be expanded in a power 
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series and one finds 

W(A lt A 2 ) - W{A X )W{A 2 ) = _ ^lMp^l + 

AMAiMr^ + n 2 2 + ni n 2 ) + A 3 (A in 2 + A 2 n 2 2 ))g\ d n x n 2 

2A 2 

6 fA 1 3 A 2 (A 2 + A 3 ) 2 (2N 3 + N)nf A 1 3 A 2 2 (A 2 + A 3 ) (2N 3 + N) n 2 n\ 
~ 92 d nm 2 ^ + + 

A X 2 A 2 2 (3A 3 (A 2 + A 3 )N 2 + A 1 (3A 3 N 2 + A 2 {AN 2 + 1))) n 2 2 n 2 

12A 3 N 

A 1 2 A 2 3 {A 1 + A 3 ) (2N 3 + N) n 2 3 m A 1 A 2 3 (A l + A 3 ) 2 (2N 3 + N) n 2 4 
+ 12A 3 N 2 + 2AA 3 N 2 

+ 0(g 7 2d ), 

(2.12) 

a result that should be reproduced by standard perturbation theory in four dimen- 
sions once we identify g 2d = —g\JA. 

The other relevant limit is of course the large N strong coupling expansion in 
which the AdS/CFT correspondence should offer the right answer. We concentrate 
our attention on the case U\ = n 2 = 1 and are interested in the normalized correlator: 
the large N limit (A = g\ d N fixed) is given as an infinite series of Bessel functions 



14, 1 



(Wi W 2 ) A 7 7 , / IA 1 A 
(W^iw^-PPA 2 ^ 2 ^ 



A X A 2 



where A, = A — Ai. In the next sections we will be interested in comparing this 
result with the prediction of supergravity. For this reason, we have to expand the 
above result for large A: the correlator in the strong coupling regime becomes 

3 



(W 1 W 2 ) A A X A 2 



(W 1 )(W 2 ) N 2 A 2 



A X A 2 IV A X A 2 



(2.13) 



The first term in the expansion corresponds to the U(l) factor present in U(N) and 
we shall drop it since it is not generally considered in the supergravity analysis. The 
first non-trivial term which can be compared with supergravity is the second one. 
This comparison is dealt with in detail in section |j. 

3. Perturbative analysis of the correlators at order g 6 

In this section we shall illustrate the main features of the numerical computation of 
the correlators of two latitudes at order g 6 (from now on we denote g^ d simply by g) . 
To be specific, we have chosen to consider two explicit configurations: 
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o symmetric CASE: The two latitudes are located at opposite positions with 
respect to the equator of the 2-sphere, namely one at 9 = 5 and the other at 
6 = n — 5, where 6 denotes the standard polar coordinate on S 2 . [See fig. [1|] 

o asymmetric CASE: The first latitude is fixed and it is chosen to be the equator 
of S 2 , while the second latitude is free to move {6 = 8 with < 5 < it). [See 

fig- a 




Figure 1: Symmetric configuration Figure 2: Asymmetric configuration 

A general remark is in order. To have the errors under control we have limited 
our numerical analysis in the range 0.7 < 5 < tt/2 for the SYMMETRIC CASE and 
for 1 < S < 7r/2 for the asymmetric case . Outside these two regions, i.e. for 
< 5 < 0.7 (symmetric case) and < 5 < 1 (asymmetric case) the requirement 
on the precision of the numerical integration becomes prohibitive for reasonable CPU 
times. 

3.1 Ladder diagrams 

To begin with, we shall consider all the diagrams which do not contain interactions. 
They can be naturally split into three families characterized by the number of field 
insertions at each latitude. Therefore, at order g 6 one has to consider the following 
possibilities 1 : g ■ g 5 , g 2 ■ g A and g 3 ■ g 3 . 

g ■ g 5 : We have four diagrams with only one propagator insertion in one of the two 
latitudes and we have schematically listed them in fig. ^. Notice that the third 
and the fourth diagram can be obtained from the first two by exchanging the two 
latitudes (C\ <-> C2) and thus we have really to compute only two diagrams. In the 
following we shall denote with t the angular parameter running over the latitude C\ 
and with s, the one spanning the second latitude C2. Then the contribution of the 

x In a diagram which does not contain interactions the power of g is simply determined by the 
number of field insertions. 
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*, A z * 3 

Figure 3: The four diagrams g ■ g 5 



diagrams in fig. |3] can be summarized as follows 

6 p 6 

S 1 • S 5 =ip P T dtl II ds i (Tr[A(t 1 )}Tr[A(s 2 )A(s 3 )A(s A )A(s 5 )A(s 6 )}) + 

JC!,C 2 i=2 

+ (Ci <- C 2 ), (3.1) 

where the symbol P in front of the integral means that the integration over the 
Si is ordered (0 < s 6 < s 5 < s 4 < s 3 < s 2 < 2n) and A stands for the usual 
effective connection constructed out of the gauge potential and the scalars. In (|3.1| ) 
the vacuum expectation value is obviously taken in the free theory and by expanding 
it in terms of free propagators we find 



i 5 $9 b N 
g g • 



4 



dti JJ dsiAi 2 (ti, s 2 ) A 22 (s 3 , s 4 )A 22 (s 5 , s 6 )+ 

JC1C2 ,_o 

(3.2) 

5 ' f 

+ -^Pf dti IJ rfs i A i2^i, s 2 ) A 22 (s 3 , s 5 ) A 22 (s 4 , s 6 ) + (Ci <-> C 2 ), 



i=2 



where Ai 2 (t;,s.,) represents a propagator connecting the latitudes C\ and C 2 , while 
An(ti,tj) and A 22 (sj,Sj) denote an internal exchange on C\ and C 2 respectively. 
Their explicit expression, if we use the polar representation for our circuits (Ci = 
{0, sin 61 sint, sin^ cost, cos^}, C 2 = {0, sin 9 2 sins, sin 9 2 coss, cos# 2 }), is given by 

^ , . sin^i sin^ ((cos^i cos^ 2 — 1) cos (tj — Sj) + sin (6*i) sin (6* 2 )) 

' 87r 2 (sin^i sin^ cos (tj — s 7 ) + cos^i cos^ — 1) /0 oN 

. 2/ , . 2o " ( 3 - 3 ) 

sin fc'i sin # 2 

Au(ti, tj) = — A 22 ( Si , Sj) = g-j-. 

The integration over the two circuits can be easily performed in a closed form for 
two generic latitudes at 9 = 9i and 9 = 9 2 and we obtain the following compact 
expression 

g 1 ■ g 5 = 9 °uAW ) {A * {A2 + Mf + MAl + A ^ A ^ ^ 
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in terms of the area A\ (A 2 ) enclosed by the circuit C\ (C 2 ) and the area A 3 delimited 
by the two latitudes. For our choice of configurations, the above expression yields 
the following two results 



g^g 5 



SYMMETRIC: 



ASYMMETRIC: 



g°(l+2jV 2 ) 
12N 



cos (f) 4 sin (|) 



5\8 



g b (l+2N 2 ) 



()j]n (ll-4cos(25) + cos(4<5))sin(§) 2 . 

g 2 • g 4 : Again we have four diagrams with two propagator insertions in one of 
the two latitudes and they are shown in the fig. f|. Using the same conventions 



(3.5) 





Figure 4: The four diagrams g 2 ■ g A 
introduced for the previous case, the contribution of the above diagrams reads 

g 2 • g 4 =1^/ dt 1 dt 2 ]Jds l (TT[A(tM(t 2 )]Tr[A(s 3 )A(sM(sM(s 6 )})o+ 

JdC 2 i=3 

+ (d <-> C 2 ). (3.6) 

The symbol P denotes, this time, both the ordering in t— integration (0 < t 2 < t\ < 
2n) and in the s— integration (0 < Sq < S5 < S4 < S3 < 2tt). If we expand the 
integrand of (|3~6| ) in terms of free propagators, we obtain 

6 



2 4 9 b N 
g g =^r- 



P(h dtidt 2 \\ dsiA 12 (ti, s 3 )A u (t 2 , s 4 )A 22 (s 5 , s 6 ) + 
i=3 



(3.7) 



+ Z/VA dt 1 dt 2 Y\_dsiA 12 (s 1 ,s 3 )A 12 (t 2 ,s 5 )A 22 (s 4 ,,s 6 ) + (C x 

JC-lC 2 i=3 

The above expression can be evaluated for generic latitudes and yields 



g 6 (N + 2N 3 ) 



(A,A 3 + A 2 A 3 + 2A X A 2 )A{A\ 



12A e N 2 
For our particular choice of the configurations this formula reduces to 

SYMMETRIC: COS (f) 2 sin (f) 10 



Co). 



(3. 



ASYMMETRIC: 



g 6 (l+2jV 2 ) 
192JV 



S\ 4 



(3.9) 



(cos (5Y - 2) sin (f ) 
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g 3 ■ g 3 : The remaining class of contributions in the absence of interaction is de- 
picted in fig H and is given by 

,6 f 3 



,3 . s 3 = g_ 

5 8 N 2 



Pcf f[ cZt i ^ + 3(Tr[^(t 1 )^(t 2 )^(t3)]Tr[^(s 4 )^(s 5 )^(s 6 )]} Q 



+ (Ct <-> C 2 ). (3.10) 
For this family of graphs it is convenient to compute separately the three different 



2r> 




Figure 5: The three diagrams g 5 • g 3 . 

contributions. The first one is similar to the diagrams considered in the previous 
cases. The sum of D 2 and D 3 in fig. |5| can be instead separated into the so-called 
abelian and maximally non-abelian part. To begin with, let us consider the diagram 
Di which is given by 

3 



9 b N 
32 



p/ f[dt i ds l+3 A 12 (t 1 ,s i )A u (t 2 ,t 3 )A 22 (s 5 ,s 6 ) + (C x ~ C 2 ). (3.11) 

J Ci Co „-_i 



'C*iC 2 i=1 

Its evaluation is straightforward and one finds 



SYMMETRIC: COS (|) 4 sin (|) 



Di 



4A 6 



(A 2 + A 3 )(A 1 + A 3 )^ 2 



ASYMMETRIC: ^ sin (<5) J sin (f) 2 . 



(3.12) 

We come now to examine the abelian part, namely the part which is separately 
symmetric in (ti,t 2 ,t 3 ) and (S4, S5,sg). We can exploit this symmetry to eliminate 
the path-ordering in the integral and to write 

3 



Ab 



9 



48N 



j> J^rft i rfs i+3 A 12 (t 1 ,S4)A 12 (t 2 ,s 5 )A 12 (t 3 ,s 6 ). (3.13) 

J Ci .-_ 1 



'C1C2 j =1 

This integral is simply the cube of the single-exchange diagram and it value is 



Ab 



6A 6 iV 2 



(3.14) 



ASYMMETRIC: sin (f) 6 . 
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Finally, we have to compute the maximally non-abelian part, whose expression is 
given by 

3 



(b [J dUds i+3 

JClC2 i= l 



A 12 (t!, s 4 )A 12 (t 2 , s 6 )A 12 (t 3 , s 5 ) + (3.15) 



+ Ai 2 (ti, s 5 )Ai 2 (t 2 , s 4 )Ai 2 (t 3 , s 6 ) + Ai 2 (ti, s 6 )Ai 2 (t 2 , s 5 )A 12 (t 3 , s A ) 
For two generic latitudes, we can perform five of the six integrations finding 



NAb 



g 6 (N 3 - N) 
: iV 2 

7rJ 



J f n da a 2 (cos 6\ cos 6* 2 — 1) (cos 6* 2 — cos 6\] 



32 7r 2 7 cos 2 Kcos^ — cos 6^) 2 + (cos^i cos6> 2 — l) 2 sin 



—-7— (cos# 2 cos^i — 1)(— 2 cos^i + cos6 l 2 (cos 0i + 2) - 1) + 2tt 3 J 3 

= g!^^) N AB[g 1 ,g 2 ], (3.16) 
where the constant J is defined by 2 

J= dsA 12 (t,s). (3.17) 
Jo 

Actually we could also perform the last integration in terms of Li 2 (z), but for the 
subsequent numerical analysis this integral representation is more useful. 
Let us collect the above results in a compact form. Apart from the maximally non- 
abelian contribution, all the other ladder graphs can be summed to give 

LAD SYM/ASYM0 = ^6 jVLAD SYM/A S YM [J] + SYM/ASYM ^ ^ (3^ 



where 



5 /A 4 /A 8 1 f6\ 2 / <P 10 



Lad5H«S] = — cos ^-J sin^-j +3Cos(^-J sin ^- j , (3.19a) 

Lad^ sym [<5] = -\- sin 2 ( - j (47- 20cos(5) - 24cos(25) + 4cos(35) + cos(45)), 
3072 V2 7 

(3.19b) 

LADj™ [<5] = cos ( ^ ) sin ( ~ ) + ^ cos ( ^ ) sin ( ^ ) + ^ sin ( ^ 
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Lad£/™[5] = sin 2 ( - J (83 - 84cos(5) + 4cos(25) + 4cos(35) + cos(4<J)). 



(3.19c) 

6144 Sm [2, 

(3.19d) 

The remaining maximally non-abelian part NAB[#i, 2 ] can be evaluated numerically 
with high precision starting from expression ( |3.16| ), with irrelevant numerical error. 

2 The result does not depend on t since A 12 (i, s) = A 12 (i — s). 
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3.2 Interaction diagrams 



We now consider all the diagrams at order g 6 containing one or more interaction 
vertices. A partial analysis of this family of graphs was performed in Jl3| and lfL5 | 
and in the following we heavily rely on the results of both papers. There it was 
shown how to reorganize the different contributions in order to get a result which is 
manifestly free of UV-divergences. In particular the diagrams were divided into three 
different classes [I Y- diagram, H- diagram and X- diagram], which are separately 
finite. 



IY-diagram This term corresponds to the sum of graphs depicted in fig. 



H they contain both the vertex contributions and the one-loop bubble corrections. 
These diagram are separately UV divergent and in order to get a finite expression it is 
convenient to collect them as illustrated in fig. R. We shall call these two quantities 




lY 



IY_ 



Figure 6: IY-DIAGRAM 



IY up and IYd own . They are finite |0| and one can obtain one from the other by 
exchanging C% with C 2 . For example the explicit expression for IY up was derived in 



and it is given by 



IY 



up 



A 3 J 
8Y 2 



2- 



dtidt 2 dt 3 ds2e(t 1 ,t 2 , t 3 ) {(ii o y 2 )2x 2 ■ d. 



y-2 



(xx o x 2 )y 2 ■ d X2 }l 1 (x 1 ,x 2 ,y 2 ) - 2 / diidt 2 ds 2 (x 1 oy 2 )l 1 (x 1 ,x 2 ,y 2 

Jo 



(3.20) 



Let us briefly recall the notation introduced in [|T^, 15|. Given two circuits x(t) 
and y(s) the effective scalar product (x o y) is a short-hand notation for x ■ y — 
\x\\y\Q± ■ Qy, where \x\Qj. = M}ei rs x r x s . Here and in the following X{ = x(ti) and 
yi = y(si) will denote points on the upper and lower latitudes respectively. The 
function I\(xi, x 2 , x%) carries the information about the integration over the position 
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of the three vertex and it is defined by 



l 1 (x 1 ,x 2 ,x 3 ) 



d 4 z 



(27r) 6 (x\ — z) 2 (x 2 — z) 2 (x 3 — z)' 



(3.21) 



One can perform the integration over z and one gets the following more useful ex- 
pression in terms of one Feynman parameter: 



l 1 (x 1 ,x 2 ,x 3 ) 



64vr 2 



da- 



1„„ ( i I ({x 3 -X2)-q(x 1 -x 2 )) 2 

lUg ^1 -1 a (l-a)(x 1 -x 2 ) 2 

((x 3 - x 2 ) - a(xi - X 2 )f 



(3.22) 



Actually one can also perform the last integration in terms of Li 2 (z), but the integral 
representation is more suitable for a numerical analysis. The expression for TY^oum is 
obtained by exchanging x with y and t with s. The final step is to evaluate explicitly 
the integration over t 3 by means of the formula dt 3 e(ti,t 2 , t 3 ) = 27rsign(t! — 1 2 ) — 
2(ti — t 2 ). If we define the function IY[#i, 9 2 ] as follows 



IY up + IY 



down 



g 6 (N 3 - N) 

W 2 



TY[0iM 



(3.23) 



its value can be computed numerically with the Montecarlo integration contained in 
Mathematica 7 both for the symmetric and for the asymmetric case. 



H-diagram I The H-diagram is drawn in fig. 0. In |L3|, [L5[ its structure was 



analyzed in great detail. For two latitudes, the contribution of this diagram can be 
cast in the following simple form 



H 



A 3 
8iV 2 



d A w P M {x 1 ,y 1 ,w)P M {x 2 ,y 2 ,w)+ 

Ai 



+ Q M (x 1 ,y 1 ,w)Q M (x 2 ,y 2 ,w) 



(3.24) 



where 



1>Z1T 

P M (x i} y h w) = dridoi [2y\ M (xi ■ d v li(x h y u w))- 
Jo 

-2xf(y i -d x / 




(2vr) 4 (xi - w) 2 (y j - wf 



Figure 7: 

(3.25) DIAGRAM 



H- 



and 



2- 



Q M (x i ,y i ,w) = Jdnda^oy^M - V) (27r )4( Xi - w )2 (y ._ w) 2 - 



(3.26) 
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In ( 3.24|) , (|3.25|) and ( 3.26| ), the index M is a ten-dimensional label running from 1 
to 10 and in particular we have defined x M = i0 7 |x|) and 8m = {9^,0). Let us 
compute first A 2 . It is convenient to rewrite this contribution as follows 

A 3 f 2lT 

A 2 = J dT 1 dT 2 da 1 da 2 x 1 oy 1 x 2 oy^(d xl -d yi ) ■ (d X2 - d y2 )H(xi, yi, x 2} y 2 ), (3.27) 
where 

1 f d 4 zd 4 w 
( i)3/i; 2)3/2) (27r) 10 J (x\ — z) 2 (yi — z) 2 (z — w) 2 (x 2 — w) 2 {y 2 — w) 2 

The action of (d Xl — d yi ) ■ (d X2 — d V2 ) on 7i(xi, y\\ x 2 , y 2 ) can then be evaluated with 
the identity (A. 7) given in j27J. One finds 

(d xi ~d yi )- (d X2 - d y2 )H(x 1 ,y 1 ;x 2 ,y 2 ) = 
1 



(x\ - yi) 2 (%2 - 3/2) 



l (4) (x 1 ,y 1 ,x 2 ,y 2 )((x 1 - x 2 ) 2 (y! - y 2 ) 2 - (a* - V2) 2 (x 2 - yi} 2 ) 

(3.29) 



+ (Y(xx, x 2 , y 2 ) - Y(yt, x 2 , y 2 ) + Y(x 2 , x%, y x ) - Y(y 2 , x 1 ,y 1 )) 



where Y{x\, x 2 , x 3 ) = l 1 (x 1 ,x 2 , x 3 )[(xi — x 3 ) 2 — (xi — x 2 ) 2 ]. Both in the symmetric 
and asymmetric case the integration over the circuits can now be carried numeri- 
cally and one determines the color-stripped contribution A2[# 1 ,# 2 ] defined by 

A 2 = 96{N3 m N) A2l9 1 ,9 2 ]. (3.30) 

Next we consider the evaluation of the A-^ contribution. This time we shall follow 
a different path in our analysis, namely we shall first perform the integration over 
the circuit analytically and then we perform numerically the integration over the 
position of the vertices. The first step is to study the function P M (w). The only 
non- vanishing components are M = 4, 5 

P\w) = - 2isin9 1 sm9 2 (cos9 1 - cos£ 2 ) {I s {9 2 )d Wo I c {9 x ) - I c (9 2 )d Wo I s (&i)) , 

(3.31a) 

P 5 (w) = - 2ism9 x sin^cos^ - cosfl 2 ) (I s (9 2 )d m I C {9 X ) - I c {9 2 )d w J s {9 1 )) . 

(3.31b) 

The function I c {5) and I S (S) are given in appendix |C]. In summary, we have to 

6 VcL 1 U. cl 1 6 

= MV-i) 1 n.Wl^.TO, (3 . 32) 

8N 2 J [w-z) 2 y ' 

Two of the eight integrations can be performed analytically (we do not present the 
cumbersome result): if we set 



- 15 - 



the remaining six integrals defining the quantity Al[#i, 9 2 ] can be computed numer- 
ically. This step is the most delicate one and the most unstable from the point of 
view of the convergence of the numerical integration. As discussed in the introduc- 
tion, the reason why we limited our analysis to the region 0.7 < 5 < tt/2 for the 
symmetric case and to the region 1 < S < it/ 2 for the asymmetric case is the 
requirement to have reliable results using the Montecarlo integration routine present 
in Mathematica 7. 



X- DIAGRAM 



There is final diagram to be considered: 
the so-called X-diagram (see fig. ||). Its expression is quite 
compact and it is given by 



X 



g 6 N(N 2 



1) 



r2\A 



2tt 



d*W 



8(4vr- r j 
'x\ o y 2 )(x2 o 



dt\dt2ds\ds2^ 

y\) - (x'i o x 2 ){y 2 ° y\) 



(3.34) 



(xi 



w 



1 {X 2 



W 



l {yi 



w) 



l (V2 



w) 



For the numerical evaluation the most convenient thing to do is 
to perform, first, the integration over the contours. Evaluating 
the integrals over the two circuits, for two generic latitudes we 
obtain the following expression in terms of I [5], I C [S] and I s [5] 
described in appendix |C| 




Figure 8: 

DIAGRAM 



X- 



x = g ^2 YN l ] j d ' w sin4 ft sin4 ^[(/(fli) 2 - Wi? - Wi) 2 ){W2?+ 

+ I s {6 2 ) 2 - I(9 2 ) 2 )} + [sin0! sin^ 2 (l - cos^ cos9 2 ){I c {6 1 )I c {6 2 ) + 

+ I s {dx)I s {9 2 )) - sin 2 B x sin 2 6 2 I(6 1 )I(6 2 )] 2 . (3.35) 

If we define the function X[6*i, 6 2 ] as 

x= 9 6 N(N*-l) x[eM (3 36) 

for our specific configurations we can proceed with the numerical integration without 
encountering particular problems. 

3.3 Comparison with YM 2 

We can now compare our numerical result with the analytic prediction given by two- 
dimensional Yang-Mills theory. Let us sum first all the contributions computed in 
the numerical analysis. The result is summarized for the symmetric case in fig. |9] 
and [TOj while the asymmetric case is given in fig. [IT] and [12]. 
The prediction in the present two cases can be derived easily from the general ex- 
pression for the correlator in the zero instanton sector given in [|T3| and reported here 
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Figure 9: symmetric case: Leading 
contribution g®N. The points are the re- 
sults of the numerical analysis, while the 
light gray line is the QCD2 prediction. 



Figure 10: symmetric CASE Sub- 
leading contribution g & /N. The points 
are the results of the numerical analysis, 
while the light gray line is the QCD2 pre- 
diction. 



Figure 11: asymmetric CASE: Lead- 
ing contribution g G N. The points are the 
results of the numerical analysis, while 
the light gray line is the QCD2 predic- 
tion. 



Figure 12: asymmetric case Sub- 
leading contribution g G /N. The points 
are the results of the numerical analysis, 
while the light gray line is the QCD2 pre- 
diction. 



in ( [2.12| ). For the symmetric case we find 



(WW) 



„, — (5 + 4cos6 + cos5 )sin f - 1 + T^y sm 1 2 ' 



(3.37) 



while for the asymmetric case we obtain 

2 



(WW) 



9 b N 



g e =—- 1 sin I — J (-36 cos(5) - 20 cos(25) + 4 cos(3<5) + cos(45) + 59) ) + 
30 72 \ V 2 



6iun 



sin - 



-52 cos(5) - 4 cos(2<5) + 4 cos(35) + cos(4<5) + 59) J . 

(3.38) 



In order to compare the results presented in figs [| [H], |TT] and [12] with the answer 
of matrix model, we compute the difference A Cal _p red between the calculated values 
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and the predicted ones. The results of this analysis are plotted in figs. |H| [11, HI 
and [Tj| where the bar denotes the estimated errors. We see that the difference from 
the central value is quite small. It is easy to see that the average absolute error is 
of order 10~ 6 / 10~ 7 , the relative error is of the order 10~ 4 at worst. 



Figure 13: symmetric case: 
Acai-Pred at the leading contribu- 
tion g 6 N 



-r+n 



Figure 14: symmetric case A Ca i_ Pred 
at the sub- leading contribution g G /N 



Figure 15: asymmetric case: 
Acai-Pred at the leading contribution 
9 & N 



Figure 16: asymmetric CASE 
Acai-Prcd at the sub-leading contribu- 
tion g 6 /N 



We also note that the error bars increase as the coincident limit is approached. This 
is a generic feature of the calculation as the integrands become increasingly singular 
in this limit. Conversely, the error is relatively small for small 5, however the preci- 
sion required to reliably calculate certain integrals in this "shrinking" limit becomes 
prohibitive for 5 less than about 0.7 (1.0), in the symmetric case (asymmetric 
case), and this defines the lower bound of our chosen range. We thus conclude that 
the conjecture is verified with a relative error of order 10 -4 in the range 0.7 < 5 < n/2 
for the symmetric case and in the range 1 < 5 < ir/2 for the asymmetric case. 



4. Correlator at strong coupling from supergravity 

At strong coupling, the AdS/CFT correspondence may be used to compute the 



- 18 - 



correlator between two Wilson loops in M = 4 SYM |28|. Practically, by taking 



the limit in which the separation of the two Wilson loops is much larger than their 
sizes, and working at large N, the correlator is computed by calculating the exchange 
of light supergravity fluctuations between the Wilson loop worldsheets. At infinite 
separation, only the lightest fluctuation modes need be considered; the subleading 
contributions stemming from the relaxation of this limit are given by the exchange 
of heavier modes. In |15[| , the contributions to the correlator from a class of modes 
(dual to the Af = 4 SYM chiral primary operators Tr($3 + i<&i) J ) which includes a 
representative of the lightest modes (i.e. for J = 2) were presented. It was found that 
the contribution of this J = 2 mode alone matched the infinite separation limit of the 
YM2 result. Beyond this, a remarkable pattern of matching between contributions 
from the modes dual to Tr($ 3 + i^^) J for general J, and corresponding terms in 
the YM 2 expression for the correlator was uncovered. This remarkable pattern of 
matching terms has since been corroborated using the techniques of localization 
fm , where it was shown that the localization conditions equate the superprotected 
operator appearing in the Wilson loop's OPE expansion discovered in |15| to precisely 
the chiral primary operator referred to above. 

Beyond this pattern of matching terms, at respectively subleading orders in the 
large separation limit, the contributions of the aforementioned dual chiral primary 
modes also include terms absent from the YM2 result. One would expect these terms 
to be removed, i.e. cancelled, by the inclusion of the other supergravity modes which 
are respectively heavier, order-by-order, compared to the dual chiral primaries. 

In fact a problem emerges before this. In a correlator calculation we are in- 
structed to sum over the exchange of all possible modes. Let us concentrate on the 
bottom of the spectrum. In addition to the mode dual to Tr($ 3 + i$4) 2 , one must 
also include the mode dual to the conjugate operator Tr($ 3 — i§i) 2 , and to the or- 
thogonal operators Tr($ 3 + $4 - $5 - $g) and Tr(3($i + $2) - l). 3 These correspond 
in the supergravity picture to various S 5 spherical harmonics of weight 2. These 
extra modes contribute at the leading order 4 in the large-separation limit, and in 
order not to spoil the agreement with YM2 should be cancelled by yet other modes. 

It happens that there are two types of further supergravity fluctuations around 
AdSs x S 5 which could potentially do the job. These are the leading fluctuations 
of the NS-NS B-field with legs in the AdS 5 and S 5 directions respectively [^|[[31] • 
They are dual to the following gauge theory operators (see appendix A of pS2"||) 

^a^b —* B-field on S 5 , 

_ (4 1) 

ip A a^ip B + 2i$ AB F+ -> B-field on AdS 5 . 

The coupling of these operators has been discussed previously in the context of the 
1/2 BPS circular Wilson loop @[fj4|- We will find that they provide leading con- 

3 Other possible modes with J = 2 do not couple to the Wilson loop, see appendix [X]. 
4 In fact the dual of Tr(3($f + $2) — 1) contributes at sub-leading order. 
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tributions of the right order and sign, but fail to cancel the offending chiral primary 
contributions due to mismatched coefficients. By going higher in the supergravity 
spectrum, we have verified that the next heaviest modes all contribute beyond the 
leading order 5 and are thus powerless to save the agreement with YM 2 . 

The interpretation of the disagreement is not clear. It could be that there is 
a problem with the supergravity limit in this instance, and that string modes are 
surviving and contributing to the correlator 6 . The strong coupling limit here could 
be subtle, since we are considering Wilson loops in the limit in which they become 
the super symmetric circles of Zarembo fl35| , where the rescaled coupling in the ma- 
trix model approaches zero f36fl . Of course, there is also the possibility that the 
YM 2 -DGRT |Tl|] Wilson loop equivalence needs to be adjusted at strong coupling, 
perhaps through the effects of the undetermined 1-loop determinant appearing in 



the localization formulae [jig] . There may also be a subtlety with the supergravity 
calculations themselves. In the remainder of this section we present these supergrav- 
ity calculations in detail, and leave the resolution of the puzzle of disagreement to 
future work. 

4.1 Preliminaries 

The fundamental string solution corresponding to the latitude DGRT Wilson loop 
was provided in jTI). We write the metric of AdS§ x S 5 as 



ds* = Nl + d^+f^±d^ + dl + cog2 ^ dnj + M2 + gin2 ^ ^ j (42) 

where the angle $ G [0, tc/2]. The worldsheet coordinates are a G [0, oo) and r G 
[0, 27r). The embedding functions are z = dfl 3 = and 

. . sin 6^ . 1 

y = smOi tanha, r = — , cp = r, x = cos#j, sini? = — -, (p — T—7T, 

cosher cosh(a + <jj) 

(4.3) 

where 9i is the position of the latitude on the sphere, i.e. its radius. Note that the 
latitude's path on the internal-space sphere is also a latitude, albeit at 

#i = \ ~ 9i (4.4) 

and so 

sin^j = — = cos 9i. (4.5) 

cosh Oi 

We would like to compute the correlator between two such latitudes at polar angles #0 
and #1, in the limit #0 — > 0, 9\ — > it. In the rest of the document we take 9\ — > iz — 9\, 
so that small 9\ indicates a latitude close to the south pole of the S 2 . 



3 There is an exception with the higher mode of the AdS^ B- field, for which a bulk-to-bulk 



propagator is not available, see section [4.5. 

6 We thank Nadav Drukker for this suggestion. 
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4.2 Dual chiral primaries 

The supergravity modes that we are interested in are fluctuations of the RR 5-form 
as well as the spacetime metric. They are by now very well known, and details can 
be found in ||§ 12| [0 |9f . The fluctuations are 



ft]' 



6 J 



7+1 ^ 



s J (x)Yj(n), 



5g a/3 = 2Jg a ps J (X)Yj(n), (4.6) 

where /i, v are AdS§ and a, j3 are S 5 indices. The symbol X indicates coordinates on 
AdS 5 and Q coordinates on the S 5 . The D^D^ represents the traceless symmetric 
double covariant derivative. The Yj(Q) are the spherical harmonics on the five- 
sphere, while s J (X) have arbitrary profile and represent a scalar field propagating 
on AdS 5 space with mass squared = J (J — 4), where J labels the representation of 
SO(6) and must be an integer greater than or equal to 2. 

The bulk-to-bulk propagator for s J is given in [28|, with normalization from fl37H . 
It is expressed in terms of a hypergeometric function 



P(X, X) = ^W J 2 F X (J, J- 3/2, 2 J - 3; 



-AW), 



W 



yu 



(4.7) 



where, 



(y ~ y) 2 + (x — x) 2 + (z — z) 2 + r 2 + r 2 — 2rr cos( 

2 3 - J N 2 J(J - 1) 



J 



B, 



Given (|4.6j ) and 
derivative, 



2TX 2 ' J Tl 2 (J+l) 2 ' 

we must construct the traceless symmetric double covariant 



(4.9) 



the details of which are given in appendix [FJ. Then, using a 10-d index M = (/i, a) 
we can express the metric fluctuations as 5qmn 



SgMN s J Yj, where 



Sg. 



II!) 



S 9 



J+l 
4 



V y 



2J (J- 1) 

^(J+T) : 



.yr 



J+l 

0<?#tf = 2 J, bg w 



d v d r + -d r 
y . 



J+l 



J+l 



r 



dl - -a, 



d; - -d y + rd r 
. V 



— r 



2J (J- 1) 

y J (J+T) : 

2 2J(J- 1) 
7 (J+l)' 



2 J sin 2 1?, 



and where we have used the fact that D 2 s J = 
the expression for the correlator 

2 

8 a X M d a X N 5g 



(4.10) 

J(J — 4)s J . We may now assemble 



(W(x)W(x)) 



(W(x))(W(x)} 




AIN 



P(X,X) 8g MN d- a X M d»X N 
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As explained at the start of this section (see also appendix |A|), at the level of 
J = 2, we have four states which couple to the Wilson loops. They correspond to 
the following scalar spherical harmonics on S 5 



^2,0 = \ cos2 sin2 ^2 e 2ilp2 , Y 2 \ Q = ~ cos 2 sin 2 2 e~ 2iip2 , 



V 2 >° ■ 



1 (3 sin 2 ■& - 1) , F 2 g = -- cos 2 cos 2tf 2 . 



(4.11) 



On the string solution we have $ 2 = 7r /2, y? 2 = 0, and so these harmonics reduce to 



^2,0 = r o 2 -2,o = Y o,o,o = 2 cos2 ^> y o 2 o°o unchanged. 

j, j 



(4.12) 



We find the following results (higher order results for Y ' ±J0 for J = 2, 3, 4 have been 



presented in [[L5| ; here we are interested in the leading order in # ; $1 which is given 
by J = 2) 



(W(ar)TV(a)) 



(W(a;))(W(x)) 
(W(a;) 



A 



i COS 2 $ 



(W(aO)(W(z)) 



8iV 2 
A 

ilw 2 



91 9\ 



0(9 



10 > 



13 age? | 3(o*0> + <%<%) | „ 



5-3 2 



2 4 



(4.13) 



where 0(9 n ) is shorthand for terms of the form 9q9\ where p + q > 10. The result 
coming from YM 2 in the large A, N limit is given by fl5fl [14] 



(W(x) W(x)) 



(W{x))(W{x)) 



A 



YM 2 



IN 2 



9 3 9 3 



0(9 



10\ 



(4.14) 



^2,2 



and so matches the contribution of one of the three modes y^+2,05 ^o 2 -2,o> ^0^0 • The 
other two modes give contributions which left uncanceled spoil the agreement with 
YM 2 . The Yq q mode contributes at subleading order, i.e. 9$9f, and so doesn't 
concern us here. In the next sections we will consider the fluctuations of the B-field 
which we will find also lead as 9q9i . However we will find that they do not remove 
the extra two contributions of the first line in (|4.13|) . 



4.3 NS-NS B-field on S 5 

Continuing up the spectrum, the next lightest modes (outside of the s J ) stem from 
the fluctuation of the NS-NS B-field which can have both legs in either the S 5 , or 
the AdS 5 directions, see eq. (2.48) and what follows it in ||29| . Here we treat the S 5 
directions, whose fluctuations correspond to an AdS 5 scalar field 



m 2 k 

a" 



4. 



(4.15) 
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The conformal dimension A of an operator related to a scalar field on AdSd+i with 
mass m is given by 



Thus here we have 



A d / , d 2 
A =2 + V m + T' 



A = k + 2. 



(4.16) 



(4.17) 



The k — 1 mode thus corresponds to a gauge theory operator of dimension 3, in the 
10 of SU(4). Consulting appendix A of P^j , we find that the operator is Sab = ipAipB- 
The antisymmetric tensor spherical harmonics Y^^(fl) obey the following equa- 
tions 

^ S %Y^ = ±2i(k + 2) (- V| B + 6) YjJJ = (* + 2) 2 Yg, (4.18) 

and may are constructed using the (regular) tensor spherical harmonics given by 

Yfep] = d a x i dpx j C [imh ... lk) x h ■■■x l \ (4.19) 

where C is antisymmetric in i, j, 1%, symmetric in I2, . . . , h, and traceless on any pair 
of indices. Using the complex basis ( A..8Q , the amoun t to a choice of sign for 

the charges associated with the angles <p, <f 2 , f3- As it turns out, our Wilson loop 

1,- 



couples only to SB^ and so we require only Y^'^. There are only two modes, given 
by 



Yj- =sin^sin^ 2 e-^ 2 , 



Y, 



1'.- 
[#,v] 



sin $ cos $2 e 



Kf.i 



(4.20) 



where we have not yet normalized the spherical harmonics. Only the first will be 
non-zero on the string worldsheet. 

The quadratic action for these fluctuations has been given in P0| , see eq. (4.3) 
therein. One has 7 



2 r , 10 , _/ 1 



2 L 



(4.21) 



ie a ^ Se B* a/3 d y B 5e ,, 



where, in units where the radius of AdS$ is unity, 1/(2k 2 ) = 4N 2 / (2tt) 5 . Subbing-in 
( |4.15|) , we find 



S — 2Ck 



4N 2 
(2vr) 5 



AdS 5 



(4.22) 



7 The leading factor of two comes because the B field is related to the A field of 1 40 by A — \/2B 
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where the constant C& encodes the normalization of the spherical harmonics. Specif- 
ically one has 

r3 



Ci 



Thus the propagator is given by 



P = ^° W A 2 Fi(A, A - 3/2, 2A - 3; -4W) 
Bk 



7T 
Y 



(4.23) 



(4.24) 



where a = (A- 1)/(2tt 2 ), and B k = 8N 2 C k /(2n) 5 . See section [12] for the definition 
of W. 

Coupling to the string worldsheet, we have 



SS 



47T 

' 2vr 



d 2 ae ab d a X M d b X N 6B_ 



MN — —1 



2tt 



(4.25) 



dadr sin $ 5B, 



where a factor of i has been included due to the Euclidean signature of the worldsheet. 
Evaluating the contribution of the k = 1 mode to the correlator we find 



(W(x) W{x)) 



(W(x))(W(x)) 



a i (elel (« + «] 



SB 



a0 



N 2 2 4 V 8 



+ o{e 8 



(4.26) 



It is straightforward to further evaluate the k = 2 contributions. They lead as 9q9\ 
and so don't concern us here. 

4.4 NS-NS B-field on AdS 5 

The supergravity action for fluctuations of the NS-NS B field with both legs in the 
AdS§ directions has been worked out in [|D|, while the dual gauge theory operator 
(for the lightest mode) has been discussed in |3T|]. The AdS/CFT correspondence 
relates linear combinations of the Ramond-Ramond 2-form potential C^ v and the 
NS-NS B field B^ v to dual operators in the gauge theory 



A = V2(B + iC), A = V2(B -iC), B 



1 



-.(A + A), C 



1 



(A- A), 
(4.27) 



2V2 2V2i 
for which the action of the modes with both legs in the AdS$ directions is given by 
S = J d 10 Xv /3^_I (V ll A vp {V tx A vp - V V A PP - V p A up ) + V Q A^V a A^) 



+ ie p p A^ v d pA T(J 



(4.28) 



8 Recall our conventions for indices: jj,, v, p, a, r, etc. denote AdS§ directions while a, /3, 7,<5, e, 
etc. denote S 5 directions. Capital roman letters denote the composite 10-dimensional index. 



- 24 - 



The equation of motion for A^ u factorizes into two first order differential equations 
(c.f. eq. (2.61) in ||), 



2k + i*D 



2(k + A)-i*D 



Ap U — 0, 



(4.29) 



where *D is the operator *DA fll> = e^ u paT d p A aT . Thus A^ v decomposes into two 
modes A\ and A 2 which obey the two first order equations respectively. In order to 
realize this at the level of the action one must introduce auxiliary fields P pu and P pv 



and write the action as []40| 

i 



S 



d 10 x, 



II [ ~-^ vpaT PiivdpA-ur + — e liVpr7T P^ u d p A aT 

IP^PV - l -V a A, u V a A^ + ie^A^dpA^ 



(4.30) 



and following another linear shift 



Ax = ^(-V a V a + 4)* A + (-V a V a + 4)"*(P - A), 
A 2 = \{~V a V a + 4)U - (-V Q V a + 4)-i (P-A), 



(4.31) 



one gets the action 



S 



d 10 x^ ( -e abcde (A lab d c A lde + A 2ah d c A 2de ) 



+ (VC-VaV* + 4) + 2)A lab A? + ( V(-V Q V Q + 4) - 2)A 2ab Af ) . (4.32) 

Expanding the fields in scalar spherical harmonics V fc , one may replace the Laplacian 
on S 5 with —k(k + 4) yielding 



V / (-V Q V a + 4) = fc + 2, A;>0, 



(4.33) 



and so A 2 is the lighter field. In fact the k = mode is not physical and can be 
gauged away (see the text underneath eq. (2.63) in |2S||). This leaves us with k — 1. 
This mode has been discussed in detail in the paper f3l| . There it is argued that the 
dual CFT operator is 

(4.34) 



B 



fj,v 



4. 4.1 Bulk-to-bulk propagator 

The bulk-to-bulk propagator for the field A 2 was given in [^TJ. The propagator is 
expressed as 

(4.35) 



V -- 



where 



2 3/2 x 



G(u) 



]ir 2 \u{u + 2)f/ 2 ' 



(4.36) 
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and where u = 1/(2W) {W being given by (|4.7|) of this document). Further, we have 

K = G', H = -(1 + u) G' - 2G, (4.37) 
prime denoting differentiation by u. The tensors T l .„ D are given by 



(4.38) 



T lu;p,u = (<^ M ) (dfiuj {dydpuj - (d u v^j (d p uj {d p d D vt^ 
d p uj (d D v)j (d u d p uj + (d u uj (^d p uj (d p d p 



= e p / x °(d p d- p u) (d x d D u) (d a 
4.4.2 Coupling to string worldsheet 



The string worldsheet couples to the B-field as per ( 4.25| ). Since our string solution 



in the AdS^ directions has only the variable (f> which depends on worldsheet-r, and 
only y and r which depend on worldsheet- a, we find 

S = / dadT ( V ' Bh + r ' B< ^ ' (4 ' 39) 
where prime denotes differentiation by o. 

We are now faced with the task of relating the fluctuations of the B-field to 
the fluctuations of the physical propagating mode A 2 . We begin by considering the 
field redefinition ( |4.31| ). The auxiliary field P a b is defined by its equation of motion 



stemming from (|4.30 ) 



hv = ^u pXa d p A Xa . (4.40) 

But, since we are interested only in the propagation of A 2 , the A field must also obey 
the first order equation of motion stemming from the first factor in ( |4.29| ), therefore 



l -e^d p A Xa = -iv (4.41) 



By (|4.31| ) we therefore have for the k — 1 mode 
V3 7 1 (i 



A 2pv = A pv - -±= (^-e pv pX(T d p A Xa - = V3A pu = VsV2B^. 



(4.42) 



The contributing k = 1 spherical harmonics are two, 

Y^ = cos sin 2 e lip2 , Y^ 10 = cos sin $ 2 e~ iip2 , (4.43) 
and each give the same contribution to the correlator 



(W(x)W{x)) 



(W(x))(W(x)) 



A ( 1 V (2tt) 5 3 

— / araa / araa cosv cosv 



(v2Vs) 



4vr 2 V x/2\/3 / AN 2 tt 3 



X 



y'y'1 :> cj>y,<t>y + r ' r ' 'P(pr;(pr + V r "P 4>y\4>f + r ' V'"P 4>r\4>y 

(4.44) 
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where we have included the factor 1/(2k 2 ) from outside the supergravity action giving 
(27r) 5 /(4iV 2 ) and the normalization of the k = 1 spherical harmonic which is 7r 3 /3. 
In the propagator ( |4.35[ ), we note that the tensor T 3 does not contribute since it 
necessarily involves a derivative by the z coordinate of ( |4.2|) , which u is independent 
of. The result evaluates to (adding a factor of two to account for the two modes in 



(W(x)W(x)) 



(W(x))(W(x)) 



r A 1 



0(9 8 ). (4.45) 



This result, in combination with (|4.26| ), does not cancel the extra two contributions 



of the first line in (|4.13|) which spoil the agreement with YM2 at the leading order. 



4.4.3 Boundary terms 

In the usual way of comparing two-point functions between supergravity and the 
CFT, the on-shell supergravity action is evaluated. However, for fields with single- 
derivative kinetic terms, like here, and also for fermions, the on-shell action vanishes 
identically. The solution has been to add boundary terms to the action. In this case 
the boundary term is BTJ |fJ2 



S 



2 1 



(4.46) 



where i,j are indices on the boundary of AdS$. The natural question arises as to 
whether the presence of such a term could affect the bulk-to-bulk correlator compu- 
tation done here. We believe it does not for the following reason. In our case the 
coupling to the boundary term is r'B^,,, but r' is zero at the boundary. Thus our 
Wilson loop has zero coupling to the boundary term. 



4.5 Heavier modes 

The modes we have considered correspond to gauge theory operators of dimension 
2 (chiral primaries) and dimension 3 (the operators ( |4.1|) ). Going one step higher in 
dimension, we have the dimension-3 chiral primaries, and at dimension-4 there are 
supergravity fluctuations of the dilaton field, massless symmetric-traceless tensor in 
AdS§ (i.e. graviton), massless AdS§ vector fluctuations (stemming from fluctuations 
of the g^ a metric components), and of course the higher KK- modes of the fluctuations 
computed here, i.e. the k = 2 modes of the B-field on S 5 and AdS$. With the 
exception of the k = 2 mode of the AdS§ B-field, where the literature provides no 
bulk-to-bulk propagator 9 , we have verified that all of these modes give contributions 
to the correlator which lead as 9fj9f. 

9 We do not expect this mode to contribute before the Q\0\ level. 
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5. Conclusions 



In this paper we have explored the relation, conjectured in between the maxi- 
mally supersymmetric M = 4 gauge theory and pure Yang-Mills theory on S 2 , in the 
zero-instanton sector. In particular, according to the localization properties of the 
four- dimensional theory established in [17], |l^], the expectation values of BPS Wil- 
son loops and their correlators should be exactly computed by some matrix model 
describing the trivial sector of the two-dimensional gauge theory. We checked accu- 
rately the conjecture at weak coupling for 1/4 and 1/8 BPS correlators of "latitude" 
Wilson loops, finding excellent agreement between Feynman diagram computations 
and the matrix model expansion at the perturbative order g 6 . At large iV and strong 
coupling we have used the AdS/CFT correspondence to test the exact expression for 
the correlator: unfortunately we were unable to find a quantitative matching with 
the matrix model expectation, even after inclusion of all the relevant supergravity 
modes. The interpretation of this disagreement is not clear and may require a better 
understanding of the strong coupling limit from the point of view of string theory 
or the subtle presence of uncanceled one-loop determinants on the field theory side. 
The resolution of this puzzle surely warrants further study. 
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A. Spherical harmonics on S 5 

We describe the metric of S 5 as follows 

ds 2 = d$ 2 + sin 2 My 2 + cos 2 $ (d$ 2 2 + sin 2 tf 2 d^\ + cos 2 ® 2 dtp 2 ) , (A.f 
where $2 £ [0, vr/2] and ip, (p 2 , <p>3 G [0, 2tt). The Laplacian is given by 



V 2 =dl - (3 tan - cot tf) d# + esc 2 ®d% 

+ sec 2 (<9 2 2 + (cot $ 2 - tan 2 ) ^ 2 + esc 2 ® 2 d% 2 + sec 2 tf 2 <9j 3 ) 



(A.2) 



The weight J scalar spherical harmonics obey V 2 Y J = — J( J + 4)Y J . This partial 
differential equation is separable and solvable. The orthogonal, but unnormalized 
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solutions are given by 
Y^ n . = w \&\ (1 + w 2 ) 1+n/2 (1 + Z 2 ) 2+J/2 e^KP+h^+jsfs) 

2F 1 (l + ±(J+\j 1 \-n), 2 + l(J+\j 1 \+ n ); l + \ 3l \, - z 2 ) (A . 3) 

2j Px(l + i(|j 2 | - |j 3 | l + i(|j 2 | + |j 3 |+n); l + |j 2 |, -t/), 

where z = tan$ and w = tan $2, and 



jie[-J,J], J -J] 1*1 = 0, 2,4, J even , J c 

n — J — Ij'iI, J - -2, |j 2 | + |j 3 |, 



J-l, J odd 
J, J even 



(A.4) 



giving the requisite (3 + J) (2 + J) 2 (l + J) /VI states, i.e. the number of components 
in a traceless symmetric rank-J tensor C^...^) in the embedding space R 6 , where the 
spherical harmonics may be expressed as 

Y J = C [h .., j) x h ...x l \ (A.5) 

where 

x 1 = sin 1? cos </?, x 2 = sin $ sin </?, x 3 = cos $ sin $2 cos<^ 2 , 

(A-6) 

x = cos $ sin $2 siny2 2 , x = cos cos t? 2 cosy2 3 , x = cos cos $2 sin</? 3 . 
The normalization of the K/f'™ ,- may be fixed using 



s 5 



jl J2J3 



2 3 (bil0 2 (b2|!) 2 r (1 + |( J - bi| - n)) r (2 + |(J - 1^1 + n)) 
* (J + 2)(n + 1) T (1 + \{J + IjiI - n)) T (2 + \{J + ^1 + n)) 



r(i + |(-|j2|- 


lis 


|+n))r(l + i(-|j 2 


1 + 


\h 


l+n)) 


r(i + |(|j 2 




73 1 


+ n)) r(l + |(|j2| 


+ 1 


73 1 


+ ra)) 



r(l + i(|j 2 |-|j3|+n)) r(l + i(|j 2 | + |j3|+n 

(A.7) 

A more convenient basis for the presentation of the scalar spherical harmonics are 
the complex variables 

zi — sintfe 1 ^, z 2 = cos ■& sin § 2 e^ 2 , -23 = cos 1? cos ~d 2 ^ ■ (A. 8) 

Using these the 6 Y 1 are given simply by {z± : z 2 , 23, z*, z 2) £ 3 }, while the 20 Y 2 may 
be summarized as 

{zf, z 2 , z%, z±z 2 , Z1Z3, z 2 z 3 , ziz 2 , Z1Z3, z 2 zl} + c.c. 

and {3|^| 2 -1, \ Z2 \ 2 -\z 3 \ 2 }. ( ' j 
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On our string solution we have $2 = tt/2 and i p 2 = y?3 = 0, which means Z3 = 0. 
However, there is a further simplification: the U(l) symmetry of the string world- 
sheets parameterized by the angle (p implies that the contribution to the correlator 
is zero unless the Y J are independent of (p. This issue has been discussed in some 
detail in |38| . This leaves the following Y 2 harmonics (normalized in accordance with 

(ED 10 ) 



*o+a = 7; cos 2 ■& sin 2 2 e 2 ** , F 2 l 2 2 = - cos 2 sin 2 2 e" 2 ** , 

\ I (A. 10) 

n 2 o°o = ^= (3 sin 2 - 1) , = -- cos 2 cos 2d 2 . 

These harmonics of the s J scalar field in (|4.6| ) correspond to the gauge theory oper- 
ators Tr($ 3 + i$ 4 ) 2 , Tr($ 3 - 2$ 4 ) 2 , Tr(3($ 2 + $ 2 ) - 1), and Tr($ 2 + $ 2 - $ 2 - $ 2 ) 
respectively. The spherical harmonics corresponding to the operators Tr($3 ± z$ 4 ) J 
for general J are 

Y*£ Jt0 = 2~ J/2 cos J sin J 2 e ±iJv2 . (A. 11) 

The 50 Y 3 are given by 

{Z1Z2Z3, zlz 2 z 3 , z\z 2 z^ z^z^z^}, 

{ziz 2 , z\z\, ZiZ 2 2 , ziz^ 2 } + cyclic permutations + c.c, 

{zl zl 4} + c.c, (A.12) 

{ Zl {\z 2 \ 2 -\z z \ 2 ), z 2 (A\z x \ 2 -l), z 3 (4ki| 2 -l)} + c.c, 

{^(N 2 + ka| 2 - 1/2), ^(2|z 3 | 2 - k 2 | 2 ), ^ 3 (2|^ 2 | 2 - N 2 )} + c.c. 

B. metric fluctuations 

The action of D^D U on a scalar field $ is, 

= - I%,dx$. (B.l) 

The Christoffel symbols for the AdS$ geometry are (comparing to (|4.2|) , here we use 
r\ — r, 0i = (j), x — r 2 cos 2 , z = r 2 sin 2 ) 

r 2 1 1 

pn _ _ r p2/ _ ' i y<t>i _ _ p</>i _ __ 





1 


p0i 


j 






1 






py 




2/2/ 



1 1 1 

p2/ _ _ pr-i py m 9 s ) 

w > L yn 1 L yy ' 1-°-^ 



where i — 1,2. The trace of D^D V $ is given by 



9 2 - 3^ + £ (V^ 2 + ^ + *. (B.3) 



10 The normalization used is f gs \Y\ 2 = 2 1 - J n 3 /((J + 1)( J + 2)). 
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C. The I functions 



wr\ 27T 

H$) = , (c i) 

a/(1 + |w| 2 - 2u> 3 cos£) 2 - A(wj + wf) sin<5 2 



US) = ^T^i f (l + H 2 -2 W 3Cos^) _ x 

4(u> 2 + u> 2 ) I ^ + _ 2u, 3 cos5) 2 - 4(w? + W D sin ^ 2 



= ^T in t f (l + H 2 -2 W 3Cos^) _ x 

4(u? 2 + u? 2 ) ^ ^(1 + H 2 - 2u> 3 cos<5) 2 - A(w\ + u> 2 ) sin5 2 

Here |u>| 2 = w\ + w\ + w| + iu|. 
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